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Abstract: In order to satisfy the higher and higher demand for communication systems, an 
Extended Binary Phase Shift Keying (EBPSK) system with very high spectra efficiency has 
been proposed. During the research, a special kind of filters was found, which can amplify 
the signal characteristics and remove utmost noise, i.e., at the point of the phase jumping 
corresponding to code "1", produce the amplitude impulse much higher than code "0", 
therefore, the aim of our study was to analyze the BER performance of the impacting filter 
assisted EBPSK-MODEM. Considering the receiver filtered "0" and 'T'signal with Rice 
amplitude distribution, just having different mean values, so the BER performance of 
EBPSK is deduced based on the classic detection theory, and compared with the traditional 
BPSK modulation both in spectra efficiency and in BER performance, which lays the 
theoretical foundation for the feasibility of Ultra Narrow Band communications based on 
EBPSK modulation. 

Keywords: Ultra Narrow Band (UNB); Extended Binary Phase Shift Keying (EBPSK); 
Rice distribution; Bit Error Rate (BER) 



1. Introduction 

With the development of the information society, there is higher and higher demand for 
communication systems. Therefore, how to transmit the information faster and better has become very 
important. Recently, high efficiency modulation technologies are receiving the attention by many 
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researchers, especially Ultra Narrow Band (UNB) [1,2] communication with its outstanding 
transmission ability in the wireless communication field. The concept of UNB communication, first 
proposed by H. R. Walker, can obtain very high data rates and high spectra efficiency, even higher 
than traditional channel capacity. Consequently, many researchers begin to pay more attention to this 
"break" performance, and the extension to Shannon's channel capacity equation has been proposed [3]. 

In wireless sensor networks [4], sensor nodes are typically powered by batteries with a limited 
lifetime, and even though energy- scavenging mechanisms can be adopted to recharge batteries through 
solar panels, piezo-electric or acoustic transducers, energy is still a limited resource and must be used 
judiciously, so efficient use of the sensor node battery's energy is an important aspect of sensor 
networks. Therefore, many researchers have paid attention to this problem and proposed many energy 
management schemes [5-7]. 

During the research, we found that as a kind of UNB system, the Extended Binary Phase Shift 
Keying (EBPSK) [8] modulation, proposed by Wu et ah, has better bandwidth efficiency [9], high data 
rates, attractive energy efficiency [10], and can greatly save energy for sensor nodes. However, the 
research on its BER performance just remains in the simulation stage and has no detailed theoretical 
derivations. So this paper will discuss this problem, which will provide the theoretical support for 
simulation results, and lay the foundation for the feasibility and rationality of UNB system. 

The rest of this paper is organized as follows. In Section 2, a scheme for EBPSK- MODEM is 
introduced, of which the detailed derivation of BER performance is given, and also the SNR 
improvement performance is analyzed in Section 3. Then, Section 4 collects some simulation results, 
and conclusions are given in Section 5. 

2. EBPSK Modulation and Demodulation 

2.1. The Definition of EBPSK Modulation 

EBPSK modulation [8] is defined as follows: 



where g 0 (?) andg^Ois modulation waveform corresponding to bit "0" and bit "1", respectively, Tis 

the bit duration, r is the phase modulation duration, and 6 is the modulating angle. Obviously, if 
r = T dXi&O = n , Eq.(l) degenerates to: 



It is just the classical Binary Phase Shift Keying (BPSK) modulation, so is named as the Extended 
BPSK. Figure 1 depicts the normalized the power spectral density (PSD) of the typical BPSK of (2) 
and the EBPSK of (1) at the same bit rate R, where bit duration T = 20T c = 20*2tt / w c , the phase 

modulation duration r for EBPSK modulation lasts one period of carrier, i.e., r = T c =2nlw c , the 
modulation angle 0 = n , the amplitude A = 1 , and the carrier frequency f c = w c 1 2n = 50KHz , 



g 0 (t) = Asmw c t,0 < t <T 




(1) 



g 0 (t) = A sin w c t, g 1 (t) = -A sin w c t, 0 < t < T 



(2) 
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Figure 1. PSD of the BPSK and EBPSK modulations, all carrier peaks are normalized 
at OdB. 




2.2. The Demodulation of EBPSK Signals 



As we know, if A 



2E h 



, the bit error ratio (BER) of the optimal BPSK receiver is calculated by: 



= Q 
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2\ 



V 



du 



(3) 



where E b is the signal energy used for transmitting one bit, and N Q is the PSD of Addition Gaussian 
White Noise (AWGN). 

However, if the same optimal BPSK receiver using traditional matched filter was utilized as the 
EBPSK demodulator, the BER performance would be much poorer since the difference in EBPSK 
modulated waveforms corresponding to "0" and "1" is very tiny and hard to detect, although in Figure 
1 more centralized PSD of the EBPSK appears. Therefore a special infinite impulse response (IIR) 
filter as given in Figure 2 is used in EBPSK receiver to produce high impulse at the phase jumping 
points r of bit Is, so as to transform phase modulation into amplitude changes [10]. And by this way it 
amplifies the signal characters as much as possible and removes utmost noise. Figure 3 depicts the 
response of this filter to EBPSK modulated signals. Obviously, a simple amplitude detector followed 
would perform the demodulation of EBPSK signals because of the existence of high impulse in coded 
1 s. Therefore, we gave this kind of filter the name of "impacting filter" [10]. 
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Figure 2. (a) The amplitude-frequency response of an "impacting filter", (b) The zoomed 
figure of (a). 
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Figure 3. The output of the impacting filter to the EBPSK modulated input. 



The transmitting EBPSK waveform 
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The output of the impulse filter 
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3. BER Performance of EBPSK Receivers 

3.1. The Derivation of a BER Expression 

In this section, the BER performance of the impacting filter based EBPSK receiver is analyzed in an 
AWGN channel. 
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The EBPSK modulated waveform corresponding to "0" is a pure sine wave, after passing the 
AWGN channel and the special impacting filter at the receiver, then the envelope r 0 of the filter output 

is with Rice distribution [11], whose probability density function (pdf) [12] is as follows: 

f „ A \ 



cr 



f r 2 + A 2 ^ 
2a 2 



r oA> 

^ 2 ) 



(4) 



where A 0 is the amplitude of the filter output, cr 2 the noise variance, and I 0 (z) the zero-order modified 
Bessel function, defined as follows: 



In 



n=0 2 



2,1 



■ n\-n\ 



(5) 



A similar analysis can also be done aiming at code "1". If we only consider those periods with 
phase jumping during z , the jumping information can be transformed into parasitic amplitude 

information after receiving filter, i.e., producing the higher amplitude impulse. At this time, its 
envelope r, still is Rice distribution, the corresponding pdf is: 



p(r l ) = — exp 



2a 2 



r A 

ycr 2 j 



(6) 



This derivation of the BER is based on the special and linear impacting filter as given in Figure 2, 
which can produce the amplitude impulse at the point of the phase jumping corresponding to code "1" 
that is much higher than the background of code "0", i.e., A, > A 0 

Let Aj = A 0 + AA = A 0 + M 0 and k > 0 . 

Therefore, assuming that code "0" and "1" be transmitted with equal probability and the decision 
threshold be U T , then the BER of the EBPSK system should be: 

= -[p(iio)+p(oii)] 
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where F 0 (x) and F { (x) is cumulative distribution function of Ricean random variable r 0 and r x 
respectively, and Q x (a,b) is Marcum's fuction, defined as follows: 



Q l (a,b) = e~ 



)/2 Z 

k=Q 



' a N 



I k (ab) 



(8) 



Sensors 2010, 10 



3829 



3.2. Calculation of the Parameters in the BER Formula 

We will now discuss in detail how to ascertain the parameters A 0 ,k ,a 2 and U T in the BER formula 

(7) in this subsection. 

(1) The amplitude Ao is evaluated as follows: 
In EBPSK modulation and via Fourier transform: 

g 0 (t) = Asmw c t < F > G 0 (w) = A—[s(w-w c )-S(w + w c )] m) 

j 

Let the frequency response of the filter beH(w) , then the signal filtered in frequency domain can 
be written as: 

n i 



Y 0 (w) = G 0 (w) -H(w) = A- [S{w -w c )-S(w + w c )] • H(w) 

j 

= A—-H(w c )- S(w -w ,)- A—- H (-w , ) • S(w + w c ) 
j J 



(10) 



So the waveform in time domain is: 

A 



y 0 (0 = F' 1 [Y 0 W] = ± [H(w c ) • - H(-w c ) • e** ] 



= —\h(w c )\ ■ e j{m) - \H(-w c )\ ■ g-^'-ft)] 
2 j 



(ID 



where q> x =ZH(w c ), q> 2 = ZH(-w c ) . The impacting filter has conjugate zero points and poles , so 
9\ ~ ~9i ' an( l |^( w c )| = |^( _vv c )| > men me received signal after filter is as follows: 

y 0 (t) = \H(w c )\-Asin(w c t + <p 1 ) (12) 

So A 0 = A ■ \H(w c )| , Aj = (l + k)- A ■ \H(w c )| , and in accordance with Eq.(3), let A = ^~jr~ ■ Finally 
we get: 



2E 

A o=\ H ^c)\-^ (13) 

(2) The variance <j 2 is evaluated as follows: 

Let the PSD of the input AWGN beN 0 /2, measured in Watt/Hz, then the variance a 2 at the output 
of the filter should be: 

ct 2 =^[\H(wfdw = N Q [\H(f)\ 2 df (14) 

(3) The threshold U T is evaluated as follows. 

Based on the binary detection model, the sketch map of the detection performance illustrated by pdf 
is shown in Figure 4, where shaded area indicates the BER. P(l I 0) and P(0 1 1) represent shaded area 
at the right and left of the threshold, respectively. Obviously, when threshold increases, P(l I 0) will 
decrease, while P(0 1 1) will increase, and vice versa. Therefore, in order to obtain a minimum BER, 
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the optimal decision threshold U T should be located at the intersection point of curve p(r Q ) and p(r { ) , 
i.e., the theoretical minimum BER should be equal to half of the area in the intersection part of 
curve p(r 0 ) and p{r x ) . 

Figure 4. The sketch map of the detection performance by pdf. 



p(r0) Threshold 




(4) The parameter k is evaluated as follows: 

According to Figure 4, the larger is the parameter k, the further is the distance of the curve 
/?(r 0 )and p(r { ) , i.e., the smaller is the area of the shaded part, i.e., the smaller is the system BER. 
Therefore, we hope to choose the larger k value in order to improve BER performance. 

Among all the parameters discussed above, how to select k , the increment in amplitude at the 
beginning of bit "1" after filtering, needs improving: 

Conceptually, the maximum of kA Q could reach the peak, say + 20dB , which is formed by a very 
narrow frequency selector and trap locating at the center of impacting filter as shown in Figure 2. In 
practice, the peak of M 0 might be limited by filter design and waveform selection (such as r , T and 
6), or even the status of channel, which still needs revealing. In fact, as an example given in Figure 3, 
where the peak of kA Q is about two times of its background, so at this situation, we may roughly get 
the estimation as k = 2 . 



3.3. Calculation ofSNRs 

In this subsection we will calculate the various signal to noise ratios (SNRs) so as to evaluate the 
EBPSK-MODEM. 

(1) SNR at the input of the impacting filter 

Similar as Ecj^14), let the PSD of the input AWGN beiV 0 /2, then from the autocorrelation 
function R(r) = ■ 8{r) we get its variance as following: 



s 2 = R(0) f ^ ■ S(0) ■ e jw 'dw - N ° 



(15) 



2n>-<° 2 2 
where the iV 0 / 2 is measured in Watt since it is the average power of the input noise, although in 
numerical value it is the same as its PSD. 
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\2E 

Let A = J - a S a i n ' th e average power of the signal at the input of the impacting filter is: 

A 2 E 

~Y = T = E b R b (16) 

which results in the SNR calculation at the input of the impacting filter as below: 

SN R in =^ = ^^ 

'" 2-s 2 NJ2 {U) 

where N Q /2 is the average power of input noise. 
(2) SNR at the output of the impacting filter 

The average SNR at the output of the impacting filter for bit "0" during T (and in theory for bit "1" 
during T -r also) is: 



E b \H(w c )\ „ E b \H{f c )\ E h R h 



SNR Q =^A 2 /a 2 = -*-• c =R 



N 



J [\H(f)\ 2 df No [\H(ffdf ^o-B neq 



(18) 



where f c is the carrier frequency, E h is the average energy for transmission of one bit, R h = l/T is 
the bit rate because both the BPSK and the EBPSK is binary modulation, so R h ■ E h is the average 
signal power transmitted, N {) I2 is still the PSD, and: 



B neq =[\H(ffdf/\H(f c )\ 2 (19) 



is the noise equivalent bandwidth of filters [13]. 

Therefore, the gain in SNR of impacting filter for bit "0" should be: 



SNR 0 _ E b R b E h -R b _ 1 



0 SNR,, N 0 • B neq I NJ2 2B neq 



(20) 



From Eq. (20) it is clear that B neq is the SNR "controlling factor" for bit "0" determined by the 
shape of frequency response of the impacting filter used. 

Similarly as Eq.(18), the peak SNR at the output of the impacting filter for bit "1" during its 
beginning is given as follows: 

A 2 A 2 -(l + k) 2 A 2 -(2k + k 2 ) 

SNR } =-^= 0 V 2 ' = 2SNR Q + 0 V 2 >- > SNR 0 (21) 

(j (j (j 

And the gain in SNR of impacting filter for bit "1" should be: 

SNR, 1 A 2 -(2k + k 2 ) 

G l = err, = ~r\ + 



SNR ; „ B u d .r\TJtf\\ 2 Af w 



U 2E b R b -[\H(ffdf 

4. Simulation 



In this section, in order to verify the correctness of theoretical analysis, we compare the BER results 
calculated by the theoretical formula with one obtained by computer simulation, as in Figure 5, where 
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system parameters in simulation are selected as: f c = 930kHz , N = 20, K = 2, A = B = 1 , and 0 = n . 
The paper chooses the receiver filter having three conjugate poles and one zero point, where the detail 
filter parameters is as ref. [10]. 

Figure 5. The BERs comparison between calculated by theoretical formula and obtained 
by computer simulation. 
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According to the above results, the theoretical result and simulation result just have 1 dB 
difference in order to obtain the same BER performance, which is caused by the detection method and 
the choosing of the decision threshold in the simulation. 

5. Discussion and Conclusions 

(1) The P e _ BPSK given in (3) is an ideal and optimum value already, while the P e _ EBPSK given in (7) 
is far from giving an optimal EBPSK receiver and still depends on the filter design. In other words, it 
leaves potential for EBPSK receivers to improve. For example, from the Figure 4 we can imagine: 

(a) With the increasing ofk, the peak distance of p(r 0 ) and p(i\) becomes large, so the shaded 
area decreases, i.e., the minimum BER decreases. 

(b) With the decrease of cr 2 , although the peak distance of p(r 0 ) and p(r v ) remains unchanged, 
their shapes are narrowing, so the shaded area, or the minimum BER, decreases also. 

Therefore, under certain conditions, the BER performance of the EBPSK modulation will 
outperform or be inferior to that of the BPSK modulation, and the spectral efficiency of the EBPSK 
system is much higher than BPSK. 
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(2) Substituting k = 2 into (21) we obtain an astonishing estimation as SN^ = ISSNR 0 , which 
implies that as long as we can obtain larger value of k by optimizing the impacting feature of the filter, 
the BER performance of the EBPSK system can be further improved. 

To sum up, aimed at a UNB system, the EBPSK-MODEM, this paper has discussed the importance 
of the special impacting filter, deduced the BER formula of EBPSK system, analyzed the SNR 
improvement performance of impacting filter, discussed the choice of optimum decision threshold, and 
analyzed the reason or possibility for EBPSK to outperform BPSK both in spectral efficiency and in 
BER performance. Research on these special filters is underway, and a detailed introduction will be 
given in forthcoming papers. 
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